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(Dated: February 1, 2008) 

Electron correlations in the two-orbital Hubbard model at half-filling are investigated by com- 
bining dynamical mean field theory with the exact diagonalization method. We systematically 
study how the interplay of the intra- and inter-band Coulomb interactions, together with the 
Hund coupling, affects the metal-insulator transition, ft is found that if the intra- and inter- 
band Coulomb interactions are nearly equal, the Fermi-liquid state is stabilized due to orbital 
fluctuations up to fairly large interactions, while the system is immediately driven to the Mott 
insulating phase away from this condition. The effects of the isotropic and anisotropic Hund 
coupling are also addressed. 



I. INTRODUCTION 

Strongly correlated electron systems with multi-orbital 
bands have attracted current interest. Typical examples 
are the manganite Lai_ x Sr x Mn03D and the ruthenate 
Si'2RuO4,0 where striking phenomena such as the colos- 
sal magnetoregistance and the triplet pairing supercon- 
ductivity were observed, stimulatinsJntensive theoreti- 
cal and experimental investigations Bo Common physics 
in the above compounds is that orbital degrees of free- 
dom together with the Hund coupling play a key role in 
realizing novel phenomena at low temperatures. 

Another interesting example is the vanadium oxide 
L1V2O4, where the heavy fermion behavior was observed 
at low temperaturesla It has been suggested that the large 
mass enhancement in this-, system may originate from 
geometrical frustrationBErcl More recently, it has been 
pointed out that degenerate orbitals are also important 
to understand the heavy fermion behavior in L1V2O4BIL2I 

It is thus desirable to discuss how the metallic 
ground state is affected by degenerate orbitals and 
also by correlations among them. Although the ef- 
fects of tlp^r^b^l^^ ge j ne p^jy-. rhave— .been explored 
extensively,BBll^e^OByS.y'BEy the stability 
of a metallic state due to orbital fluctuations has not been 
studied systematically. In this paper, we study the two- 
orbital Hubbard model at half-filling to discuss how the 
interplay of the intra-band interaction, the inter-band in- 
teraction and the Hund coupling affects the stability of a 
metallic phase against an insulating phase. In particular, 
by using dynamical mean- field theory (DMFT), we show 
that orbital fluctuations stabilize a heavy fermion state 
when the intra- and inter-band Coulomb interactions are 
nearly equal. 

This paper is organized as follows. We introduce the 
model Hamiltonian and briefly summarize the formula- 
tion based on DMFT in Sec. O In Sec. 151 we then dis- 



with the exact diagonalization method, and clarify how 
the intra- and inter-band couplings affect the stability of 
the Fermi liquid phase. The effects of the Hund coupling 
are addressed in Sec. IV. A brief summary is given in 



cuss the metal-insulator transition by combining DMFT 



Sec. 

II. MODEL AND METHOD 

A Two-orbital Hubbard Hamiltonian 

We consider a correlated electron system with two-fold 
degenerate orbitals, which is described by the following 
Hubbard Hamiltonian, 

H ^ ^ ^ij^iaa^ aa ^ ^ ] ^iotl^ial 

<i,j> ia 

+ U' / ] nii^naa' +J}Sj\- S i2 , (1) 

1,(7,17' 1 

where c\ a(7 (ci a(7 ) creates (annihilates) an electron with 
spin <t(=T, I) and orbital index a{= 1, 2) at the ith site, 
and riiaa = c\ aa Ci aa . The corresponding spin operator 

is defined by S ja = \ Yjoo' c ]a<7 T <™' c j<x<r'' where r is the 
Pauli matrix, iy represents the transfer integral, U and 
U' the intra-band and inter-band Coulomb interaction, 
and J the Hund coupling. 

In the following, we will particularly focus on how the 
interplay of U, U' and J affects the metal-insulator tran- 
sition. 



B Dynamical mean-field theory 

We make use of DMFT,E1E1IEE1 which has 
been extensively .. ri.|ed for the single-band ..Hubbapi 
modeimfflffl the two-band modeLBMiSEl 
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the periodic Anderson model,@0'@ etc. This method 
has also been applied to the degenerate Hubbard, 
model by combining it with the exact ,diaaonalization,E3 
quantum Monte Caria simulation,E2rLirc30 and iterative 
perturbation theoryBJ 

In DMFT, the lattice model is mapped to an effec- 
tive impurity model, where local electron correlations are 
taken into account precisely. The lattice Green function 
is then obtained via self-consistent conditions imposed 
on the impurity problem. The treatment is exact in 
d — » oo dimensions, and even in three dimensions, DMFT 
has successfully explained interesting physics such as the 
Mott metal-insulator transition. We introduce here the 
impurity Anderson Hamiltonian with two orbitals, 

-Hhnp = J]] tkc\ aa Ckuo + Effgcrfata 
k,a,a aa 

+ J]] Vfe r<4 ao ./ acr + faa c ka<rj 
k : a.cr 

+ U ^2 f a ffaTf al fal + U' ^ flaflvfla'hv' 



JSi • S2, 



(2) 



where /1 CT (/ QCT ) creates (annihilates) an electron with 
spin a in the a-th orbital at the impurity site, and 
S « = lE TO '/i<rW/^' Note that the effective pa- 
rameters in the impurity model, such as the spectrum 
of host electrons et and the hybridization Vk , should be 
determined self-consistently so that the obtained results 
properly reproduce the original lattice problem. This 
will be done explicitly below. We focus on the symmet- 
ric case with half-filled bands by setting Ef — — \U—U', 
for simplicity. 

To discuss the stability of the Fermi liquid state in a 
normal metallic phase, we introduce the Green function 

with two components ^ = f^j as 



G(r 



-( 9 \T*(t)&(t')\ 9 ), 



(3) 



where T is the time-ordering operator and \g) represents 
the ground state. In the non-interacting case U = U' = 
J = 0, the Green function, Gq = GqI, is given by 



Go 1 (*) = z 



T v k 2 

^ z-e k 



(4) 



When the interactions are introduced, the full Green 
function G& (z) is written as, 



Z — Sintra(^) — £fc 

— Si n tor(^) Z ' 



~ ^intcr(^) 
Sjntra(^) — Cfc 



(5) 



We have here used the fact that the self-energy 53 is inde- 
pendent of the momentum in d — > 00 dimensions. Then 
the diagonal element of the local Green function is given 
as, 



Gloc (z) 



= r — [*-w*) -*]*(*) ^ (6) 

J— 00 (z — Eintra(-2j — x ) ~ Winter (2) 

where p{x) = J2k $ ( £ fc — x ) ^ s the density of states in 
the non-interacting case. The lattice structure is not 
so important to discuss the metal-insulator transition 
in d — 00 dimensions. We use here the semicircular 
density of states p(x) = -^y/l — (x/D) 2 , which corre- 
sponds to an infinite-coordination Bethc lattice. From 
the Dyson equation G^ 1 = G 1 + 53, we then obtain the 
self-consistent equation as, 



Go 1 (z) 



Gloc (z) ■ 



(7) 



We wish to note that the self-energy does not appear 
in these self-consistent equations, allowing us to simplify 
the iteration procedure. Namely, by estimating the diag- 
onal element of the local Green function for the Anderson 
model eq. (||) and by using the self-consistent eqs. (j|) 
and (0), we can discuss Fermi-liquid properties in the 
system. In the following, we take the band width D as 
the unit of the energy. 



III. METAL-INSULATOR TRANSITION 

To solve the self-consistent equations in DMFT, we use 
the exact diagonalization method. The diagonal element 
of the Green function for the Anderson Hamiltonian eq. 
(H) is given as, 



G(z) = (g\f 



1 



n 



z + Eq-H 
t 1 



+ {9lf ^ z-E + H 
f(z + E ,f^\g)) 
-f(-z + Eo,fn\g)). 



fh\g) 
At Iff) 



(8) 



where Eq is the ground state energy. The function 
/ (uj, is given by the following continued fractionCa, 



LO — tl 



LO — 0,1 



(9) 



lo - a 2 



where a n (b n ) is the nth diagonal (subdiagonal) element 
of the Hamiltonian H tridiagonalized by the Lanczos 
method with the initial vector \ip). To solve the self- 
consistent eqs. (|j) and (]?]) in terms of the obtained Green 
function G, we introduce the following cost function \ 2 
proposed by Caffarel and KrauthEa as 



1 ^ 

n=0 



|[GS ld (<«„)] 1 -[GT'(i Wn )r 1 (10) 
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where n max is chosen sufficiently large, and u> r , 



-f-*> 

is the Matsubara frequency. In terms of the conjugate 
gradient method, we determine the parameters tk and 
Vk in Gq cw by minimizing \ 2 m ec l- (P^t)- To properly 
converge the iteration in DMFT, we introduce the "ficti- 
tious" inverse temperature /?, which is fixed as (3 = 50 in 
this paper. The number of site is set as N = 6. We note 
that a careful scaling analysis for the fictitious tempera- 
ture and the number of sites is necessary only when the 
system is near the metal-insulator transition point. To 
discuss the stability of the Fermi liquid state, we define 
the wave-function renormalization factor in terms of the 
self-energy around u> = 0, 



Z = 



1 - 



AlmS (iui) 



Auj 



(11) 



which corresponds to the weight of a quasiparticle exci- 
tation. 



A Interplay of intra- and inter-band interactions 

Let us first discuss the Fermi liquid properties of a 
metallic state in the absence of the Hund coupling (J = 
0), and clarify how the orbital-degeneracy affects a metal- 
insulator transition. By making use of the exact diago- 
nalization method, we iterate the procedure mentioned 
above to obtain the results within desired accuracy. 

We show the quasi-particle weight Z as a function 
of the interband interaction U' in Fig. |l|. In the case 




FIG. 1: The quasi-particle weight Z as a function of 
the inter-band Coulomb interaction U', obtained by the 
exact diagonalization (N — 6) within DMFT. 

U = 0, Z decreases monotonically with increasing U', 
and a metal-insulator transition occurs around U' c ~ 0.9. 
On the other hand, when U ^ 0, there appears non- 
monotonic behavior in Z , namely, it once increases on the 
introduction of the interband Coulomb interaction, has 
the maximum value in the vicinity of V ~ U, and finally 
leads to a metal-insulator transition at a critical value 
of U' . It is easily understood that the large value of U' 



suppresses hopping among sites, thereby causing a sharp 
decrease in the renormalization factor Z. However, it 
is somehow unexpected that the maximum structure ap- 
pears around U ~ U' , and furthermore is more enhanced 
for larger U and U'. We think that this is related to or- 
bital fluctuations, which may stabilize a normal metallic 
state around U ~ U' ' . 

In order to observe the above behavior more clearly, we 
show the contour plot of the quasi-particle weight in Fig. 
0. When W = 0, the system is reduced to the single-band 




U' 



FIG. 2: The contour plot of the quasi-particle weight Z 
in the case J = 0. The bold-dashed line represents the 
phase boundary of the metal-insulator transition, which 
is deduced by estimating the values of U and U' that 
give Z = 0. The solid circles are the transition points 
obtained by the linearized DMFT. 

Hubbard model, where the intra-band Coulomb interac- 
tion U induces the metal- insulator transition at U c , which 
was estimated, by the numerical renormalization groiia 
(U c = 2.94) ,E3 the exact diagoiiaiization (U c = 2.93) ,E3 
the linearized DMFT (U c = 3).l3 These values agree with 
the results in Fig. || well. 

There are several remarkable features in this phase di- 
agram. We first notice that the value of Z is not so sen- 
sitive to U' for a given U (> U'), except for the region 
U ~ U'. In particular, the phase boundary indicated by 
the dashed line in the upper side of the figure is almost 
flat for the small U region. An important point to be 
noted is that when U ~ U' the metallic state is remark- 
ably stable against a transition to the insulator, and per- 
sists up to fairly large Coulomb interactions. Moreover, it 
becomes immediately unstable, once the parameters arc 
away from the line U — U' . This tendency becomes more 
conspicuous in the regime of strong correlations. These 
remarkable properties about the metallic phase around 
the line U = U' are closely related to orbital fluctua- 
tions, as we will see momentarily in the following. 

The second point is that there are two insulating 
phases corresponding to the regions U > U' and U < V , 
which are separated by the metallic phase for small U, U', 
but are continuously connected to each other via the re- 
gion of large U, U'. For example, let us observe phase 
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transitions with U = 4.0 being fixed. When U' = 0.0 
with U = 4.0, the system is in the insulating phase caused 
by the intra-band interaction U . Introducing the inter- 
band coupling U' , a phase transition occurs to the metal- 
lic phase, and further increase in U' induces the second 
transition to another insulating phase, which is domi- 
nated by the inter-band interaction U' . These two in- 
sulating phases should show different properties, though 
they can be adiabatically connected to each other in prin- 
ciple. 

B Critical properties around the phase boundary 

To clarify the above characteristic properties around 
the metal-insulator transition, we exploit a linearized ver- 
sion of DMFT,c3 which provides us with a transparent 
view of the phase transitions. In this scheme, we intro- 
duce a specific Anderson impurity model connected to 
only one host site, which corresponds to the model eq. 
(^J) with k = 1. This extreme simplification even provides 
sensible results as far as low-energy properties near the 
transition point are concerned. Namely, for low-energy 
excitations around the Fermi surface, the Green function 
for /-electrons may be approximated by a single pole as, 

GiocCO ~ p (12) 

where the residue w corresponds to the quasi-particle 
weight. Recall that for w = the ground state belongs 
to an insulating phase, while for m / to a metallic 
phase. We can thus determine the metal-insulator tran- 
sition point from the condition w = 0, as has been done 
before. By combining the self-consistent eqs. (Q) and 

(a) (b) (c) 

■5.++ + * 

E f ++ — * 

1 2 

FIG. 3: Three distinct singlet states in the Anderson 
model, which is hybridized with only one effective host 
site. The numbers 1, 2 specify the orbital indices. 

, we iterate the linearized DMFT. The self-consistent 
condition imposed on the hybridization now reads 

V=^. (13) 

For the single-band Hubbard model (U 1 — 0), the critical 
value lX-,= 3 was already estimated by the linearized 
DMFT.EI 

In this simplified model, the metallic ground state is in 
the spin singlet state drawn in Fig. p| (a) schematically. 



Although the introduction of the interband Coulomb in- 
teraction U' lowers the states (b) and (c) energetically, 
these states do not affect the ground-state properties up 
to a certain critical value of U'. In fact, the relevant ex- 
citation gap is independent of U' in this approximation, 
and the residue of the Green function around V = is 
given as, 

u> = 36QQ , (14) 

via a straightforward calculation. Then the self- 
consistent eq. dl3] ) updates the hybridization via 
(3/U)V — > V in each iteration process. Therefore, when 
U > U c (= 3), the effective hybridization V vanishes by 
the iteration, thereby stabilizing the insulating phase in 
this region. On the other hand, the hybridization is rel- 
evant for U < U c and thus the metallic phase is realized. 
Note that the critical value t/ c = 3 obtained here is the 
same as that for the single-band Hubbard model.CJ We 
thus come to the conclusion that the inter-band Coulomb 
interaction has little effect on the metal-insulator tran- 
sition, which is essentially described by the mechanism 
for the single-band Hubbard model, as far as the case of 
U > U' is concerned. This statement is indeed consistent 
with the numerical results shown in Fig. ^, namely, the 
phase boundary is rather insensitive to the change in U' 
for small U'. 

On the other hand, in the special case U = U', the 
nature of the transition is totally changed in this approx- 
imation. In this case, the ground state is formed by three 
degenerate states shown in Fig. |^, and not only spin but 
also orbital degrees of freedom play an important role. In 
fact, quantum fluctuations among these states decrease 
the effective mass, making the metallic phase more sta- 
ble. The critical value U c = 5 estimated by the present 
scheme agrees with the numerical one in Fig. § Also 
this is consistent with the conclusion of quantum Monte 
Carlo simulations JH3 which claims that the ground state 
in the single (double) orbital Hubbard model belongs to 
the insulating (metallic) phase in the case U = U' = 4. 

The present analysis also gives a clear picture for two 
types of the metal-insulator transitions: the insulating 
phase in the upper region of in Fig. ^ corresponds to the 
state of Fig. |^(a) , where the intra-atomic repulsion gives 
rise to the insulating phase, whereas the insulator in the 
lower region, corresponding to Fig. |](b) an d (c), is sta- 
bilized by the inter-band interaction. Around the critical 
point on the U = U' line, the metallic state sandwiched 
by these insulating phases is realized by the competition 
among almost degenerate singlet states (a), (b) and (c). 
In other words, the metallic state is stabilized by orbital 
fluctuations. 



IV. EFFECTS OF THE HUND COUPLING 
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A Isotropic case 

We now discuss the effects of the exchange coupling J 
between orbitals, which plays an important role in real 
materials such as manganitcs and ruthcnatcs. By per- 
forming similar calculations, we obtain the contour plot 
of the quasi-particle weight Z as shown in Fig. We 
also plot the values of Z as a function of J in Fig. H. 




FIG. 4: The contour plot of the quasi-particle weight 
Z in the system with the ferromagnetic Hund coupling 
between the orbitals J = —1.0. 



Z0.6 




FIG. 5: The quasi-particle weight Z as a function of J. 

The Hund coupling J (< 0) enhances spin correlations 
between the orbitals, thus renormalizing electrons near 
the Fermi surface. It is indeed seen in both figures that 
for U > U' the quasi-particle weight is decreased with 
the increase of |J|. On the other hand, the effect of the 
Hund coupling is less conspicuous in the region U < U'. 
In the latter region, the inter-band Coulomb interaction 
dominates electron correlations in the ground state when 
J = 0, which is schematically drawn as Fig. || (b) or (c). 
For these configurations the Hund coupling is irrelevant, 
which explains that the effect of J is less important for 
U < U'. 
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FIG. 6: The density of states for V = 0.0, 1.2, 1.6 and 
2.0 obtained by the exact diagonalization with S = 0.01 



To make the metal-insulator transition clear, we also 
calculate the density of states by means of the exact di- 
agonalization. The results are shown in Fig. ||. When 
U = 1.0, J = —1.0 and U' = 0.0, the system belongs to 
the metallic phase, as shown in Fig. |]. Introducing the 
interband coupling V, fluctuations between the orbitals 
are enhanced, which assist to stabilize the quasi-particle 
states up to U' ~ 1.2. Further increase in U' favors the 
doubly occupied state in the same orbital. Therefore, the 
electron states around the Fermi surface are pushed away 
to higher energy regions, as seen in the case of U' = 1.6, 
and finally a quantum phase transition occurs from the 
metallic phase to the insulating phase. In the insulating 
phase for U 1 = 2.0, we can clearly see the charge gap in 
the density of states in Fig. ^. 

We also investigate the case of the antiferromagnetic 
exchange coupling (J > 0). The contour plot in the case 
J = 0.7 is shown in Fig. u\. By comparing with the ferro- 
magnetic case, it is seen that electrons are renormalized 
strongly, as shown in Fig. ^. This difference is explained 
by counting the number of effective degrees of freedom at 
each site in the presence of the exchange coupling. In the 
ferromagnetic case, the system has a tendency to form 
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FIG. 7: The contour plot of the quasi-particle weight Z in 
the system with the antiferromagnetic exchange coupling 
between the orbitals J = 0.7. 



the S = 1 state at each site, while in the antiferromag- 
netic case it favors the singlet state made of two orbitals. 
The resulting effective degrees of freedom give rise to the 
difference in the stability of the Fermi liquid state: the 
metallic state with the ferromagnetic coupling is more 
stable than the antiferromagnetic case. Note that the 
situation is similar to the stability of the metallic state 
due to orbital fluctuations around U ~ U' discussed in 
the previous section. 



B Anisotropic case 

We finally discuss the effects of anisotropy in the Hund 
coupling. For this purpose, we introduce the anisotropy 
parameter A as, 

(S 1 -S 2 ) A = A(5^|+5f5|) 

where A = corresponds to the Ising-t 
while A = 1 the isotropic case. In Fig 
quasi-particle weight Z for U — 1.0 and J = —1.0. In 
the ferromagnetic case, the quasi-particle weight Z is lit- 
tle affected by anisotropy although somewhat nonmono- 
tonic dependence is found as a function of A as shown in 
the inset. On the contrary, the anisotropy has a notice- 
able effect for the antiferromagnetic case with small U', 
as seen from Fig. |[ where the anisotropy decreases the 
weight w monotonically with the increase of A. This dif- 
ference arises from the fact that the x — y components of 
the exchange coupling bring about quantum fluctuations 
more prominently in the antiferromagnetic case than in 
the ferromagnetic case, giving rise to the sizable effect on 
the antiferromagnetic case in the small U' region. 



SfSI, (15) 

e anisotropy, 
we show the 




FIG. 8: The quasi-particle weight Z for the case with 
anisotropic exchange coupling: U — 1.0 and J = —1.0. 
Inset shows the quasi-particle weight as a function of the 
anisotropy parameter A. 




FIG. 9: The quasi-particle weight Z as a function of U' 
in the case U — 1.0. 



V. SUMMARY 

We have investigated the stability of a metallic phase 
in the two-orbital Hubbard model at half-filling by means 
of dynamical mean-field theory with particular emphasis 
on the role played by orbital degrees of freedom. By com- 
bining the exact diagonalization method with DMFT, we 
have discussed how the interplay among the intra-band 
and inter-band Coulomb interaction together with the 
Hund coupling affects the metal-insulator transition. In 
particular, it has been found that the metallic state is 
remarkably stable, even up to U ~ 5.0 if the intra and 
inter-band Coulomb interactions are nearly equal. Also, 
slight deviation from this condition immediately drives 
the system to the insulating phase. We have pointed out 
that orbital fluctuations play a particular role to realize 
the metallic state around U ~ U' ~ 5.0. The effect of the 
exchange coupling has been also discussed. It has been 
clarified that the effective degrees of freedom at each site 
play an important role again in stabilizing the metallic 
state. 
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We have focused on nonmagnetic phases of the model 
in this paper. Possible instabilities toward other ordered 
states such as the magnetic order, the superconductiv- 
ity, etc., are to be studied. In particular, the ferro- 
magnetic instability is interesting in connection with fer- 
romagnetism realized in manganite compounds. These 
problems are now under consideration. 



VI. ACKNOWLEDGMENTS 

We would like to thank M. Suminokura for useful dis- 
cussions. This work was partly supported by a Grant- 
in- Aid from the Ministry of Education, Science, Sports 
and Culture of Japan. A part of computations was done 
at the Supercomputer Center at the Institute for Solid 
State Physics, University of Tokyo and Yukawa Institute 
Computer Facility. 



REFERENCES 

1 Y. Tokura, A. Urushibara, Y. Moritomo, T. Arima, 
A. Asamitsu, G. Kido and N. Furukawa, J. Phys. Soc. 
Jpn. 63, 3931 (1994). 

2 Y. Maeno, H. Hashimoto, K. Toshida, S. Nishizaki, T. 
Fujita, J. G. Bcdnorz, F. Lichtenberg, Nature, 372, 
532 (1994). 

3 M. Imada, A. Fujimori and Y. Tokura, Rev. Mod. 
Phys. 70, 1039 (1998); Y. Tokura and N. Nagaosa, 
Science 288, 462 (2000). 

4 T. M. Rice and M. Sigrist, J. Phys. Condens. Matt. 7, 
L643 (1995). 

5 S. Kondo, D. C. Johnston, C. A. Swenson, F. Borsa, 

A. V. Mahajan, L. L. Miller, T. Gu, A. I. Goldman, M. 

B. Maple, D. A. Gajcwski, E. J. Freeman, N. R. Dilley, 
R. P. Dickey, J. Mcrrin, K. Kojima, G. M. Luke, Y. 
J. Uemura, O. Chmaissem and J. D. Jorgensen, Phys. 
Rev. Lett. 78, 3729 (1997). 

6 H. Kaps, N. Buttgen, W. Trinkl, A. Loidl, M. Klemm 
and S. Horn, J. Phys. Condense. Matter, 13, 8497 
(2001). 

7 M. Isoda and S. Mori, J. Phys. Soc. Jpn. 69, 1509 
(2000). 

8 S. Fujimoto, Phys. Rev. B 64, 085102 (2001). 

9 H. Tsunetsugu, preprint. 

10 Y. Yamashita and K. Ueda, preprint. 

11 H. Hasegawa, J. Phys. Soc. Jpn. 66, 1391 (1997). 

12 S. Ishihara, M. Yamanaka and N. Nagaosa, Phys. Rev. 
B 56, 686 (1997); R. Maezono, S. Ishihara and N. 
Nagaosa, Phys. Rev. B 58, 11583 (1998). 

13 R. Fresard and G. Kotliar, Phys. Rev. B 56,12909 

(1997) . 

14 A. Klejnberg and J. Spalek, Phys. Rev. B 57, 12041 

(1998) . 

15 J. Biinemann, W. Weber and F. Gebhard, Phys. Rev. 
B 57, 6896 (1998). 



16 Y. Motome and M. Imada, J. Phys. Soc. Jpn. 67, 3199 
(1998). 

17 G. Khaliullin and S. Maekawa, Phys. Rev. Lett. 85, 
3950 (2000). 

18 T. Takimoto, Phys. Rev. B 62, R14641 (2000). 

19 T. Momoi and K. Kubo, Phys. Rev. B 58, R567 
(1998). 

20 M. J. Rozenberg, Phys. Rev. B 55, R4855 (1997). 

21 K. Held and D. Vollhardt, Eur. Phys. J. B 5, 473 

(1998) . 

22 J. E. Han, M. Jarrcll and D. L. Cox, Phys. Rev. B 58, 
R4199 (1998). 

23 V. S. Oudovenko and G. Kotliar, Phys. Rev. B 65, 
075102 (2002). 

24 G. Kotliar and H. Kajuetcr, Phys. Rev. B 54, R14221 
(1996). 

25 W. Metzner and D. Vollhardt, Phys. Rev. Lett. 62, 
324 (1989). 

26 E. Muller-Hartmann, Z. Phys. B 74, 507 (1989). 

27 A. Georges, G. Kotliar, W. Krauth and M. J. Rozen- 
berg, Rev. Mod. Phys. 68, 13 (1996). 

28 T. Pruschke, M. Jarrell, and J. K. Freericks, Adv. 
Phys. 42, 187 (1995). 

29 M. Caffarel and W. Krauth, Phys. Rev. Lett. 72, 1545 
(1994). 

30 Th. Pruschke, D. L. Cox and M. Jarrell: Phys. Rev. 
B 47 (1993) 3553. 

31 O. Sakai and Y. Kuramoto, Solid State Comm. 89, 
307 (1994). 

32 R. Chitra and G. Kotliar, Phys. Rev. Lett. 83, 2386 

(1999) . 

33 J. Joo and V. Oudovenko, Phys. Rev. B 64, 193102 
(2001). 

34 M. S. Laad, L. Craco and E. Miiller-Hartmann, Phys. 
Rev. B 64, 195114 (2001). 

35 R. Bulla, Phys. Rev. Lett. 83, 136 (1999). 

36 A. Georges, G. Kotliar and W. Krauth, Z. Phys. B 
92, 313 (1993). 

37 Th. Maier, M. B. Zolfl, Th. Pruschke and J. Keller, 
Eur. Phys. J. B 7, 377 (1999). 

38 Y. Imai and N. Kawakami, J. Phys. Soc. Jpn. 70, 2365 
(2001). 

39 Y. Ono, R. Bulla and A. C. Hcwson, Eur. Phys. J. B 
19, 375 (2001); Y. Ohashi and Y. Ono, J. Phys. Soc. 
Jpn. 70, 2989 (2001). 

40 M. Jarrell, H. Akhlaghpour and T. Pruschke, Phys. 
Rev. Lett. 70, 1670 (1993). 

41 T. Mutou and D. Hirashima, J. Phys. Soc. Jpn. 63, 
4475 (1994); 

42 T. Saso and M. Itoh, Phys. Rev. B 53, 6877 (1996); 

43 R. Bulla and M. Potthoff, Eur. Phys. J. B 13, 257 

(2000) . 

44 R. Haydock, V. Heine and M. J. Kelly, J. Phys. C 8, 
2591 (1975). 



